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Abstract Quantum field theoretic treatments of fermion oscillations are typically restricted to calculations 
in Fock space. In this letter we extend the oscillation formulae to include more general quasi-free states, and 
also consider the case when the mixing is not unitary. 

1. Introduction 

The theoretical underpinning for fermion oscillations was developed decades ago [6], but only recently 
' have terrestrial and solar neutrino experiments begun to substantiate this work [22,23], and experimental 

evidence for oscillations has led to renewed interest in the theory. Current quantum field theoretic treatments 
of the fermion oscillation phenomena [1,2,7-18,21] have brought about modifications to the oscillation formula 
developed in [6], adding terms depending on the sum of energies rather than their difference. Although these 
new results are based on calculations made in fermion Fock space, they rely on modifications suggested by 
physical considerations. In this paper we generalise the oscillation formula to general quasifree states, and 
show that the additional terms occur naturally in that setting. Our result contains the known formula for 
• Fock states as a special case, but also includes other physical scenarios such as the thermal (KMS) state, or 

situations in which polarisation of the vaccum has occurred. The oscillation formula for the thermal state 
could be a better approximation for fermions at nonzero temperature such as solar neutrinos. At the end of 
the paper we shall also consider the case when the mixing operator is not unitary and show that it leads to 
. similar effects. 

| From a mathematical point of view the main obstacle to such calculations lies in the fact that we wish to 

. calculate the oscillatory behaviour of correlations between projections onto flavour states (such as the v e and 

states) at different times, but the dynamical behaviour is simplest in states with definite masses. These are 
distinct from the states of definite flavour to which they are related by a non-trivial mixing transformation. 
The fermionic anticommutation relations (CAR) can be written in terms of smeared creation operators 

c(w) = J w(x)a*(x)dx 

(smeared with test functions w, z in the complex inner product space of wave functions on R 3 with values 
in the product of Dirac spinors V and an iV-dimensional space V describing the various flavour states) as 

•i— I , 

[c(w)*,c(z)] + = (w,z), [c(w),c(z)]+ = 0. 

5-H 



It is well-known [19] that the (abstract) algebra defined by these relations has many inequivalent representa- 
tions by operators in Hilbert space, and the study of the interrelations between a selection of these forms the 
main focus of this paper. We shall mainly be concerned with quasi-free representations of the CAR algebra 
which generalise the standard Fock and Dirac-Fock representations. (The Fock representations themselves 
have been studied, for example, in [7,8,20,24].) Quasi-free representations are those in which the Wick de- 
terminant formula expresses the n-point correlation functions in terms of the 2-point correlation functions 
just as in Fock space. As well as appearing for thermal states of systems, quasi-free states often arise in 
situations where the vacuum is polarised, and so allow us to treat more complicated field-theoretic effects 
whilst avoiding the detailed models. A well-known technique of Powers and St0rmer [25] and Araki [3] tells 
us how to construct any quasi-free representation of CAR(PF) as the composition of an injection of W into 
W © W with a Fock representation of CAR(IF © W), and so we shall concentrate on Fock representations. 

2. The one particle space 

To establish notation we first recall that for Dirac particles the elements of W = W m can be thought of 
as the initial data for the Dirac equation 

in-— = H D w 
at 
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where the Dirac Hamiltonian Hp is given in terms of momentum operators P by 

H D = c(a ® P + (3 ® Mc) 

with a = (ai, a 2 , 0^3) and /3 satisfying the Clifford algebra relations 

aj/3 + (3<Xj = 0, 1 — 1, Oj-afc + CKfeCKj = 2<5jfc, j, fc = 1, 2, 3, 

and M a positive operator on V with eigenvalues the masses mi,f/i2, . . . , ra„. (In what follows, we shall 
omit the tensor produt and write 0M for /3® M, etc.) 

Choosing a basis in which the mass matrix M is diagonal, we refer to the solutions of the Dirac eqaution 
as mass eigenstates. We define 

P± = i (1 ± HdE- 1 ) , 

where E is the positive square root of the positive operator Hp = (\P\ 2 +M 2 c 4 ). (We shall also use Ej for the 
value when in the eigenstate with mass rrij.) The P± are idempotent, self-adjoint, and P±Hr> = EP±; that 
is, they are the positive and negative energy projections on W. They also determine the mass representation 
of the CAR algebra with creation operators c m and a Dirac-Fock vacuum vector O m which satisfies the Dirac 
condition that, for every w in W, 

c m (p + w)*n m = = c m (p_w)n m . 

When we study the flavour space Wf we work rather with a flavour representation cy in a standard 
Fock space with a flavour vacuum f2/ satisfying Cf(w)*Clf = for all w e Wf. 

3. General mixing transformations 

In our previous paper we studied what happened when the mass and flavour spaces W m and Wf were 
isomorphic by a unitary mixing transformation T, but here we shall consider more general situations such 
as orthogonal mixing transformations, and Powers-St0rmer transformations which enable us to realise quasi- 
free states on Fock spaces. When T is unitary the spaces W m and Wf can be identified, but, for orthogonal 
T, when W m and Wf are the same as real spaces but have different complex structure, it is simpler to treat 
them as distinct. 

It will be convenient to consider more generally the case when we have two inner product spaces Wj 
and Wk and a map Tjk : Wk Wj which is orthogonal in the sense that it preserves the real part of the 
inner product: for all w, z £ Wk, 

(z,w) + (w,z) = (T lk z,Tj k w) + (Tj k w,T jk z). 

For any real-linear operator Tjk on Wk we define the complex linear map ajk — \{Tjk — JjTjkJk) and the 
antilinear map bj k = j{Tjk + JjTjkJk) where Jj and Jk simply indicate multiplication by i on Wj and Wk, 
respectively. The Fock space creation and annihilation operators Cj and Ck are linked by 

c k (w) = Cj(aj k w) + Cjibjkw)*, 

and 

Cj(T jk w) + Cj(T jk w)* = Cj(a jk w) + Cj(b jk w)* + Cj(a jk w)* + Cj(b jk w), 

for all w E Wk, where ajk and bjk are the Bogoliubov maps just defined. To be consistent with the 
anticommutation relations in both Wj and Wk we have the orthogonality relations 

a *jk a ik + b* jk bj k = 1 = a lk a* k + b jk b* jk , a* jk b jk + b* k a jk = = a jk b* k + b jk a* k . 

These are just the conditions that T be orthogonal. Since these transformations mix creation and annihilation 
operators it is expedient to introduce a more succinct notation. We combine creators and annihilators in 
the row vector Cj = (cj c* ) , and introduce 

A ( a 0k b ]k \ 

3k ~W a jk )> 
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to obtain 



c k {w) = Cj(A jk w), for w 
One advantage of dealing with the general situation is that there is an obvious composition law 

Aji = AjfeAfej, 

obtained by writing the relationship between creation operators on Wj and Wi directly and through the 
intermediate space W k - The orthogonality properties of the Bogoliubov maps can also be interpreted as 
telling us that a jk = a* k] , and b jk = b* k] , or 

(This is closely related to Araki's self-dual construction [3].) 

Instead of the vacuum states described above we shall use a more general quasi-free state cu for the CAR 
algebra. As noted above, this is determined by the two-point correlation functions which define a complex 
linear operator R and a conjugate linear operator S by 

ui[c(w)* c(z)] — (w,Rz) 
w[c(w)c(z)} — (Sw,z) 

where R = R*, < R < 1, and S = -S*, [3,4,5]. We note that if K = i(2R - 2S - 1) defines a complex 
structure on W, that is K 2 = — 1, then the state w is a Fock state for some choice of complex structure. 

The GNS construction guarantees the existence of a representation 7r : CAR(14 7 ) — > TL containing a 
cyclic vector VL € Ji such that cu(b) = (fi, 7r(6)0) for any b e CAR(W). The representation tt over H can 
be expressed in terms of a Dirac-Fock representation of CAR(W + where W + and are both 

isomorphic to W. The representation in question takes (w + , W-) G W + ® W- to c + (w + ) + c_(w_) and the 
u> + with a Dirac type vacuum which is killed by the annihilators c + (w + )*, and the creators c_(w_). This 
means that 

(c + (w)Q,c + (z)Q) = (w,z), (c_(w)*^,c_(z)*0) = (z,w), 

and all other two-point correlation functions vanish. Suppressing the representation map, and writing c = 
iroc, the required representation of CAR(VF) is given by 

c(w) = (c+(a + w) + c + (b + w)* + c_(a_w) + c_(&_u>)*) 

where a± are linear and b± are antilinear and satisfy the orthogonality relations a±a± + b±b± = 1, a±b± + 
b*±a± = 0. (It is easy to check that this does provide a representation of CAR(I4 / ).) 
When W — Wj we shall write a±j and b±j for the Bogoliubov maps to give 

C J'( W ) = 72 ( C +( a +J W ) + C +( h +3 W Y + C -( a ~j W ) + c -i h -J w T) ■ 

We adopt the summation convention that when an index s can range over different values {+,— } and is 
repeated (as in expressions such as c s (a s jiv) or Cj(aj S a s kw)) one sums over all its values and divides by y/2. 
Then the above expression can be abbreviated to Cj(w) = c s (A s jw), and it is easy to see that our earlier 
rules for compositions apply. Thus a state which gives a quasi-free representation for the label j will also do 
so for the label k. 

With this notation we have Cj(w) = c s (A S jw). When applied to 0, only certain of the components are 
non-zero. For instance, 

Cj(w)fl = -j= {c+(a +j w) +c+(b +j w)* +c-{a-jw) +c-(b-jw)*)Cl= -j= (c+(a +j w) + c-(b-jw)*)n, 
and, similarly, 

Cj (w)*Sl = -j= (c + (b +3 w) + C-(a- jW )*) 0. 
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Reversing the expansion we obtain 

Cj(w)fl = i (cj((a j+ a +j + bj-b-j)w) + c j ((a j _6_ j + b j+ a +j )w)*) 0, 

and 

Cj(u;)*n = - (cj((a j+ b +j + bj-a-j)w) + Cj((aj-a-j + b j+ b +j )w)*)fl, 



which may be combined as 



Cj(w)il = Cj(A?w)tt, where A° = ^ S R/ ^ 



where R — ^{a 3+ a + j + bj-b-j), S— ^(aj+b+j + 6j_a_j), 5* = i;(aj-b-j +bj + a + j) (which the orthogonality 
relations show is the adjoint of S), and R' = ^(dj-d-j + bj + b + j) = 1 — R (by the orthogonality relations). 
In effect A° is the projection on Wj associated with Q. We easily calculate that 

(c j (z)Q,c j (w)n) = (cj(z)n,(cj(Rw) + Cj(Sw)*)Q) = {z,Rw) 

so that R is the correlation operator already introduced, and similarly for S. 

Whenever U is a unitary transformation of Wj, there is a unitary map itj(U) implementing U in the 
sense that 

niWciiwfriiU)- 1 = cj(Uw), 

for all w € Wj. Since iTj{U) is unitary there is a similar relation for the annihilator Cj(w)* and we combine 
these as 

niWcjiwfriiU)- 1 =cj(Uw), 

where; 

U=( U ° 

\o u 

(When j refers to flavour this gives a representation of the group of flavour transformations.) 

We shall also need the infinitesimal version of this which arises from taking the unitary exp(isP) and 
differentiating at s = to obtain 



[n'j(P),Cj(w)} = -i- ^{e isP )cj{w)^j{e isP )-\ =Q = -i- Cj(e isP w)\ s=Q = c 3 (Pw). 



When P = 1, we can regard N = tt'j(1) as the number operator, and when P = P x , the projection onto the 
states of flavour A, then N\ = tt'^(P x ) counts the number of flavour A particles. 

For our calculations we need to know the effect of iTj{U) on Ck(w) for different j and k. We therefore 
note that 

TTjiUJCkiw^jiU)- 1 = TTjiU^jiAjk^TTjiU)- 1 

= c 3 (UA jk w) 
= c k (A kj UA jk w) 

= c k ({akjUajk + b k jUb 3k )w) + c fe ((b k jUa jk + a kj Ub 3k )w)* . 
It is convenient to introduce the abbreviation U k = A k jUAj k , and write 

■Kj(U)c k (w)-Kj(U)~ x =c k (U k w), 
(it being understood that U is an operator on W k ). We have the explicit formula 



U k = 



u k v k 
v k u k 



with 



Uk — akjUajk + bkjUbjk, ffe — akjUbjk + bkjUdjk- 



4. Correlation functions 

The correlation functions which we first wish to calculate have the form 

(c f (z)n,ir m (U)*Tr' f (P)ir m (U)c f (w)n), 

where P projects onto a flavour state, U — exp(—itHjj/?i) gives the time evolution, and f2 defines a quasi-free 
state. This can easily be found by differentiating the more tractable 

{cf(z)Sl,ir m (U)*Wf{D)ir m {U)cf(w)Sl), 

where D = exp(isP). 

Gathering together our various comments we calculate that 

ir m (U)*ir f (D)Tr m (U)c f (w)n = n m (U)*7r f (D)n m (U)c f (A f w)n 

= Z f (U* f DU f A f w)TT m (U)*n f {D)Tr m {U)n. 

We now differentiate this to get 

*7r m (I7)*7^(P)7r m (l7)c/(ii;)n = c f (U}iPU f A° f w)n + icf(U}UfA f w)n m (U)*Tr' f (P)Tr m (U)fl 

= c f (A° f U}iPUfA° f w)n + ic f (A f ti)ir m (U)*n' f (P)ir m (U)n. 

Inserting this expression into the inner product, but with w = ^ ™ j , we obtain 

(c f (z)n, ir m (U)*ir' f (P)ir m (U)c f {w)Sl) = -i(c f (z)n,c f (A f U}iPU f A° f w)n) 

+ (c f (z)n,c f (A° f w)n m (Uy7r' f (P)7r m (U)n). 

In each inner product we take the creation operators from the right to an adjoint acting on the left. There 
the factor of A° ensures that the adjoint annihilates ft, so that we simply get 

(cf{z)Q,n m {Uy^ f (P)n m {U)c f (w)n) = -i{[c f (A ^* f iPUfA f w)*,c f (z)} + n,n) 

+ ([c f (A f w)*,c f (z)} + n,7r m (U)*7T' f (P)7r m (U)n). 

The anticommutators can now be written explicitly in terms of inner products. For example, the second 
gives 

[c f (A° f wr,c f (z)} + = [c f (Rwy+c f (Sw),c f (z)} + = (Rw,z). 
The first requires a more detailed calculation, but we require only the first entry in A^U^iPU f A®w: 

<> »)(? t)C« £)(:; z)(s %)(" 



<* **>(;? :|)(o p ° P ){7, :;)(?:) 



Recalling that v is conjugate linear, the product of the three middle matrices can be written as 

u* f v}\ ( iP 0\(u f v f \_.(F G 
v* f u* f ) \ iP) [v f Uf) - l \G F 



where 

F = u)Pu f - v}Pv f , G = u* f Pv f - v}Pu f . 

Using the fact that R = R* and recalling the conjugate linearity of S, this enables us to rewrite the first 
commutator as 

(R S*)i(^g p^^g™^=i(R*FR + R*GS-S*GR-S*FS)w. 

Combining the expressions for the two commutators we obtain 

(z, (R*FR + R*GS - S*GR - S*FS)w) + (z, Rw){Sl, 7r m (C/)*7^(P)7r m (f7)n). 

The physically interesting quantity is the expected number of flavour A particles (Nx(t))^, in a state 
where one flavour fi particle has been created out a time t earlier out of the "vacuum" f2. This can be 
obtained by taking P = P x , z = w = P^(f>j, U — V t — exp(—iHDt/K) above, and then summing as <f>j runs 
over an orthonormal basis of W, to get (with F x and G x denoting F and G when we take P — P x ) 

J2(cf(P<*4> 3 )n, ir m {Vt)W f (P)n m {Vt)c f {Pi l 4> J )Sl) = tr [(R*F X R + R*G X S - S*G X R - S*F X S) P"] 

3 

+ tr [RP»] (n,7T m (V t )*TT' f (P x )7r m (V t )n). 

The last inner product is just the vacuum expectation of P x and should be subtracted (since we are only 
interested in the enhancement produced by creating a flavour state), and we must also divide by the norm 
of the state tr[RP»] to get 

_ tr [(R*F X R + R*G X S - S*G X R - S*F X S) P^] 
" K[RPT] ' 

The expected total flavour number (N(t))^ is obtained by summing over A, (which means that F x and 
G x are replaced by u*jUf — vjvf and u^Vf — VfUf, respectively. The ratio (N\(t)) M / (N(t)) M then gives the 
proportion of flavour A particles. In the next two sections we shall look at two special cases of this formula. 

During the preparation of this paper an interesting preprint appeared [17] which investigates the CP 
violation in three flavour mixing. We note that in our context the transition probability for antiparticles 
can be calculated by using annihilation in place of creation operators, which leads to replacement of R by 
R' = 1 - R and of S by S* = -S, to give 

_ tr [((1 - R)*F X (1 - R) - (1 - R)*G X S + S*G X (1 - R) - S*F X S) P^] 



tr [R'Pv] 



5. Unitary mixing in quasi- free states 

The first case which we shall consider is for a unitary mixing matrix but in a quasi-free state. When the 
mixing transformation is given by a unitary operator T, we have Uf = TUT* and Vf — 0, so that G x = 
and F x = TV t *T* P x TV t T* . The expression then simplifies to 

tr \(R*TV t *T*P x TV t T*R- S*TV t *T*P x TV t T*S) P»] 
= • 

(This formula can easily be checked in the case of a thermal state ojp at absolute temperature (fc/3) -1 , where 

k is Boltzmann's constant (see the Appendix), and gives the known values R = (l + e^ 1311 ®) and 5 = 0, 
in agreement with our general formula.) 
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We note that summation over A gives 

W)), 



tr [{R*TV t *T*TV t T*R - S*TV t *T*TV t T* S) P»] _ tr[(R*R - S*S)P^} 



tr[i?P^] tr[RP»] 



which is, as one would hope, independent of time 

In fact for states invariant under global U(l) transformations of W, we always have S = [4,5], and 

then 

, . ti[R*TV t *T*P x TV t T*RP^} 

{Nx(t)) « = ^RP^} • 

To achieve a more explicit formula we take R = Tp(H]j)T* to be a function of the Hamiltonian, where 
p(x) is a real function of the real variable x, defined everywhere except, perhaps, zero. Both the Dirac and 
KMS states satisfy this restriction. This gives 

, tr[Tp(H D )V*T*P x TV t p{H D )T*P»] tr[p{H D )V*T* P x TV t p(H D )T* P»T] 

( ~~ tr[TpT*Pn] ~ tr[pT*PKT} 

Performing some preliminary calculations, we have 

p(h d )v; = p(H D y H ^ h (p + + p_) 

= p{E)e lEt ' n P+ + p{-E)e- lEt ' n P- 

p(E)e iEt/h + p(-E)e- tEt/h ] + \ \p(E)e iEt/n - p(-E) e - lEt/h ] (P+ - P_), 



and, in particular, 

T*RT = \ [p(E) + p(-E)] + | [p(E) - p(-E)} (P + -P.). 
To condense our notation, we define 

g j = p{E i )e iE i t l h + pi-E^e-W* 
5j = p{Ej)e iE ^ t/n - p{-Ej)e- iE i t/h 
jj = p{Ej) + p{-Ej) 
ej = c(a ■ P + (3mjc)/Ej. 

We note that ej has trace zero but 

tr(cje fc ) = (|P| 2 c 2 + mjmkC^/iEjEk) = S jk . 
We also know that in terms of the mass basis the components of P x are 

(T*P x T) 3k = T* x T xk = T Xj T xk . 
With this notation the numerator in the oscillation formula is 

AT N 

\ ^[(<J 3 + S 3 e 3 )T*P x T(a k + S k e k )*T*P^T] = Wj*k + 5jS k S jk ]T Xj T xk T ^ k T M 

j^k—l 

N 

= \ ^ [(fj^fc + ^"^fe)(l + Sjk) + {ojOk - Sj5 k ){l - S 'j k )]T \jT Xk T 'nkTpj . 
j,k=i 

Now, recalling that p is a real function, we have 

aj<j k + Sjh = 2[p{E 3 )p{E k )e^ E ^ E ^' h + p(-E j )p(-E k )e i ^- E ^ h ] 
a 3 a k - 6jS k = 2[p{E 3 )p(-E k y^+ E ^l h + p(-E j )p(E k )e- i ^ +E ^ h ]. 
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When 771 j = TTik we have Sjk = 1 so that the second term in the numerator disappears, and since also 
Ej = Ek, we see that <Jja k + SjSk is time- independent, and consequently there is no flavour oscillation 
between these. For Fock states and general masses, one has p(E) = 1 when E > 0, and p(E) = when 
E < 0, which gives UjO k - SjS k = 0, and UjU k + SjS k = 2e l{E ^ Ek)t/n . 
In general, the denominator is 

N 

tr [r] = 2^ 7j W, 

i=i 

Whenever p(E) + p(—E) = 1, as happens for Fock states and also the thermal states where p(E) = (1 + 
e~ f3E )~ 1 , we have 7j = 1, and then the unitarity of T means that the denominator is 2, giving 

N 

{N\(t))n = \ ^ [( a i w k + 8j8 k )(l + Sj k ) + {ojOk - SjS k )(l - Sjk)]T\jTx k T IJl kTfj l j 
j,k=i 

and our earlier formulae for ajlik ± SjSk show that this contains both standard oscillations depending on 
the energy differences and others depending on Ej + E k ■ In the Fock vacuum state the oscillation formula 
is consistent with the calculations performed in [20]. 

Taking A = p, the oscillation formula takes on a slightly more compact form 

N 

(Nn(t))n = 5 [( a j°k + SjS k )(l + S jk ) + {djtJk - 5j5 k ){l - Sjk)] |T W -| 2 |T Mfe | 2 . 
j,k=i 



6. Non-unitary mixing in a Fock state 

We could instead work with the flavour vacuum. Then there is no need to inject W into W © W, so 
that we have S = and R = 1 This gives 

tr \F X PH 
(N x (t))^ = — L__i, 

where 

F x = u}P x u f -v* f P x v f , 

with 

Uf = a fm V t a m f + bf m V t b m f, vj = af m V t b m f + bf m V t a m f. 

From this we may show that both sorts of oscillation terms occur in this case too. However, the total flavour 
number is given by replacing F x by F — u*^uj ~ VjVf, and 

<*<*». = 

and in general this depends on the time t. This is essentially a squeezing phenomenon. It provides a strong 
reason to be cautious about non-unitary mixing of this kind. 

Appendix 

For a thermal state at temperature (fc/3) -1 the KMS condition and anticommutation relations give 
formally 

ujp[c(w)*Dc(z)} = ujp [Dc(z)c(e^ HD w)*] 

= ujp[b {{e (iIiD w, z) - c{e fmD wYc{z))] 

= ( e P HD w,z)uj p [D}-LUp[c(eP HD w)*Dc(z)}+LUp[c(DeP HD w)*c(z)}. 



This can be rearanged as 

Uf3{c((l + eP HD )w)*Dc(z)} = {e fSHD w,z)uj [D]+uj [c{De pHD w)*c{z)], 
or, replacing w by (1 + e l3HD )^ 1 w, 

u [c{w)*Dc(z)] = <(1 + er^^w, z)u p [D] + u fj [c{D(l + e- pBD )- l w)*c(z)]. 
The case D = 1 (and D = 0) gives the usual two point correlation function 

up[c{w)*c(z)] = ((lle-^rV*), 

so that 

W(3 [cH*5c(z)] = ((1 + e-^)- 1 ^, «)w /J [5] + <(1 + e- 0HD )- l D{l + e~ pHD )- l w, z). 



Acknowledgements. We are grateful to M. Blasone and Tsou S.T. for useful conversations on different 
aspects of neutrino oscillation. 

The second author would like to thank the Rhodes Trust for their support. 

References 

E. Alfinito, M. Blasone, A. Iorio, and G. Vitiello, 'Squeezed Neutrino Oscillations in Quantum Field 



[1] 

[2] 

[3 

M 
[5] 

[G] 

[7] 

[8] 

[9] 

[10 
[11 

[12 
[13 

[14 

[15 

[16 

[17 



Theory 'Phys. Lett. B, [hep-ph/9510213 . 

E. Alfinito, M. Blasone, A. Iorio and G. Vitiello, 'Neutrino Mixing and Oscillations in Quantum Field 



Theory', Acta Phys.Polon. B 27 (1996) 1493-1502, hep-ph/9601354 



H. Araki, 'On quasifree states of CAR and Bogoliubov automorphisms', Publ. RIMS Kyoto Univ., 6 
(1970-71), 385-442. 

E. Balslev and A. Verbeure, 'States on Clifford algebras', Commun. Math. Phys. 7 (1968), 55-76. 

E. Balslev, J. Manuceau and A. Verbeure, 'Representations fo anticommutation relations and Bogoliubov 

transformations', Commun. Math. Phys. 8 (1968), 315-326. 

S.M. Bilenky and B. Pontecorvo, 'Lepton mixing and neutrino oscillations', Phys. Rep., 41 (1978), 
225-261. 

M. Blasone and G. Vitiello, 'Quantum field theory of fermion mixing', Ann. of Phys., 244 (1995), 
283-311. 

M. Blasone and G. Vitiello, 'Erratum: Quantum field theory of fermion mixing', Ann. of Phys., 249 
(1995), 363-364. 

M. Blasone, P.A. Hen ning and G. Vit iello, 'Mixing Transformations in Quantum Field Theory and 
Neutrino Oscillations', hep-ph/9605335 . 



M. Blasone, 'New Results in the Physics of Neutrino Oscillations', hep-ph/981032£ . 
M. Blasone, P.A. He nning, and G. V itiello, 'The exact formula for neutrino oscillations', Phys. Lett. B 
451 (1999) 140-145, |hep-th/9803157 . 



M. Blasone, P.A. Henning, and G. Vitiello, 'Green's Functions for Neutrino Mixing', b.cp-ph/9807370 



M. Blas one and G. Viti ello, 'Remarks on the neutrino oscillation formula', Phys. Rev. D 60 (1999) 
111302, |hep-ph/9907382 . 

M. Blasone, P. Jizba, and G. Vitiello, 'Currents and charges for mixed fields', Phys. Lett. B 517 (2001) 
471-475, |hep-th/01036"87 . 

M. Blasone, A. Capolupo, and G. Vitiello, 'Comment on "Remarks on flavor-neutrino propagators and 
oscillation formulae" ', bep-ph/0107183[ 

M. Blasone, A. Capolupo, and G. Vitiello, 'Understanding flavor mixing in Quantum Field Theory', 
hep-th/0107125| . 

M. Blasone, A Capolupo and G. Vitiello, 'Quantum field theory of three flavour neutrino mixing and 
oscillations with CP violation', |hcp-th/0204184. 



9 



[18] K. Fujii, C. Habe, and T. Yabuki, 'Remarks on flavor-neutrino propagators and oscillation formulae' 



Phys.Rev. D 64 (2001) 013011, hep-ph/0102001 



[19] L. Garding and A.S. Wightman, 'Representations of the commutation and anticommutation relations', 

Proc. Nat. Acad. Set. (1954) 40, 617-626. 
[20] K.C. Hannabuss and D.C. Latimer, 'The quantum field theory of fermion mixing', J.Phys A: Math. 

Gen., 33 (2000), 1369-1373. 
[21] C-R. Ji and Yu. Mishchenko 'The general theory of quantum field mixing', Phys. Rev. D 65 (2002), 

096015-40. 

[22] E. Kearns, T. Kajita and Y. Totsuka, 'Detecting massive neutrinos', Scientific American, 281 (1999), 
48-55. 

[23] S.F. King, 'Neutrino oscillations: Status, prospects, and opportunities at a neutrino factory', J. Phys. G: 

Nucl.Part. Phys., 27 (2001), 2149-2170. 
[24] R. Plymen and P. Robinson, Spinors in Hilbert space, Cambridge University Press, Cambridge, 1994. 
[25] R.T. Powers & E. St0rmer, 'Free states of the canonical anticommutation relations,' Commun. Math. 

Phys. 16 (1970), 1-33. 



10 



